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Current  density  distributions  and  local  state  of  charge  (SoC)  differences  that  are  caused  by  temperature 
gradients  inside  actively  cooled  Li-ion  battery  cells  are  discussed  and  quantified.  As  an  example,  a  cylin¬ 
drical  Li-ion  cell  with  LiFeP04  as  cathode  material  (LiFeP04-cell)  is  analyzed  in  detail  both  experimentally 
and  by  means  of  spatial  electro-thermal  co-simulations.  The  reason  for  current  density  inhomogeneities 
is  found  to  be  the  local  electrochemical  impedance  varying  with  temperature  in  different  regions  of  the 
jelly  roll.  For  the  investigated  cell,  high  power  cycling  and  the  resulting  temperature  gradient  additionally 
cause  SoC-gradients  inside  the  jelly  roll.  The  local  SoCs  inside  one  cell  diverge  firstly  because  of  asymmet¬ 
ric  current  density  distributions  during  charge  and  discharge  inside  the  cell  and  secondly  because  of  the 
temperature  dependence  of  the  local  open  circuit  potential.  Even  after  long  relaxation  periods,  the  SoC 
distribution  in  cycled  LiFeP04-cells  remains  inhomogeneous  across  the  jelly  roll  as  a  result  of  hysteresis 
in  the  open  circuit  voltage.  The  occurring  thermal  electrical  inhomogeneities  are  expected  to  influence 
local  aging  differences  and  thus,  global  cell  aging.  Additionally  the  occurrence  of  inhomogeneous  current 
flow  and  SoC-development  inside  non-uniformly  cooled  battery  packs  of  parallel  connected  LiFeP04-cells 
is  measured  and  discussed. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

According  to  literature  the  most  important  impact  factors 
on  calendar  life  aging  of  Li-ion-cells  are  temperature  and  SoC 
[1-3],  the  main  influencing  factors  on  cycle  life  are  temperature, 
current  density,  depth  of  discharge  per  cycle  (ADoD),  and  SoC 
[1,3-5].  Generally,  the  electrochemical  aging  effects  lead  to  an 
increase  in  cell  impedance  and  a  loss  of  cell  capacity.  Recently, 
Liu  et  al.  [6]  reported  that  the  aging  progress  of  LiFeP04/graphite 
systems  that  are  also  investigated  in  this  work  is  dominated  by 
the  loss  of  active  lithium  and,  thus,  induces  mainly  the  loss  of 
capacity  with  only  little  changes  in  the  cell  impedance.  The  data 
of  capacity  loss  presented  in  that  article  show  also  a  dependency 
on  temperature  and  current  rate  for  cycled  LiFeP04/graphite 
systems.  Most  cycle  life  models  (e.g.  [4,5])  consider  global  values 
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of  one  or  more  cell  states  and  cycling  conditions,  such  as  tem¬ 
perature,  current  density  and  SoC  [1-5].  This  approach  neglects 
thermal  and  electric  inhomogeneities  inside  the  cell;  the  impact 
of  temperature  gradients  and  the  resulting  current  density  and 
SoC  gradients  across  the  cell’s  jelly  roll  is  not  considered.  The 
methods  and  results  presented  in  this  article  contribute  to  a  better 
understanding  of  macroscopic  electric  and  thermal  interactions 
during  cycling,  which  are  influenced  by  the  cell  design  and  the 
cooling  design  outside  the  cell.  It  is  beyond  the  scope  of  this  article 
to  parameterize  a  spatialized  aging  model  for  LiFeP04  batteries. 
However,  to  reproduce  the  effects  that  mainly  influence  the  aging 
processes  (namely  current,  state  of  charge,  and  temperature)  a 
thermal-electric  cell  model  with  spatial  resolution  is  introduced 
in  this  work.  Thus  temperature  gradients,  current  density  distribu¬ 
tions,  and  gradients  of  SoC  inside  a  battery  cell  during  cycling  can 
be  simulated.  This  detailed  cell  simulation  is  of  high  interest  for 
interpreting  electric  and  aging  behavior  of  Li-ion  cells,  which  are 
not  only  influenced  by  the  maximum  temperature  of  the  cell  but 
also  by  the  temperature  distribution  occurring  inside  electrode 
stacks.  This  novel  approach  provides  a  valuable  tool  for  examining 
the  temperature  spread  inside  battery  cells  and  modules  for  system 
design. 
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Fig.  1.  Pulse  impedance  ZP(10s)  of  the  investigated  LiFePC^-cell  after  10s  of  high 
current  pulses  for  discharge  and  charge  loads  (88  A). 
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Fig.  2.  Ratio  of  10  s-discharge  pulse  impedance  to  10  s-charge  pulse  impedance  for 
the  investigated  LiFePC^-cell  dependent  on  temperature  and  SoC. 


2.  Theoretical  background 

For  a  better  understanding  of  cell  state  inhomogeneities,  first 
the  macroscopic  thermal-electric  effects  that  are  expected  to 
occur  inside  the  jelly  roll  of  a  Li-ion  cell  during  operation  should 
be  highlighted.  In  this  section  the  theoretical  background  is 
described  qualitatively  without  using  simulative  results.  A  vali¬ 
dated  thermal-electric  simulation  model  reproducing  these  effects 
is  introduced  in  Section  3.  All  figures  in  this  section  that  contain  val¬ 
ues  of  electric  cell  performance  are  based  on  measurements  on  the 
cylindrical  LiFeP04  cell  used  in  the  experimental  section.  Its  electric 
and  thermal  characteristics  can  be  seen  in  Table  1. 

2.1.  Temperature  gradient 

During  cycling  Li-ion  cells  generate  heat  mainly  in  consequence 
of  overpotential  under  current  load  (irreversible  heat)  and  changes 
in  entropy  (reversible  heat).  This  thermal  energy  heats  up  the 
cell  and,  as  cell  temperature  exceeds  ambient  temperature,  it  is 
conducted  to  the  cell  surface,  where  it  is  dissipated  into  the  envi¬ 
ronment.  As  a  result  a  gradient  in  temperature  from  the  warmer 
inner  elements  of  the  jelly  roll  to  the  colder  outer  elements  occurs 
inside  the  cell  [7].  The  simulative  description  of  heat  generation  and 
of  the  occurring  temperature  distribution  inside  the  cell  is  declared 
in  Section  3. 


2.2.  Gradient  of  current  density  and  heat  generation 


In  order  to  compare  the  electric  behavior  of  a  Li-ion  cell  for  dif¬ 
ferent  temperatures  in  the  time  domain  the  term  pulse  impedance 
Zp(t)  is  used.  It  describes  the  cell’s  normalized  overvoltage  response 
on  constant  current  steps  in  dependence  on  time. 


W  = 


A  U(t) 
A  I(t) 


(1) 


Fig.  1  shows  the  measured  values  ofZp  after  10  s  pulse  load  for 
the  investigated  cell  as  a  function  of  temperature.  All  values  depend 
on  measurements  at  a  starting  SoC  of  50%. 

Considering  the  thermal-electric  interdependencies,  we  can 
assume  that  warmer  regions  of  a  Li-ion  cell  comprise  lower  elec¬ 
tric  impedance  than  colder  cell  regions  and— as  the  regions  are 
electrically  connected  in  parallel  to  the  colder  cell  regions— exhibit 
higher  current  densities  during  operation.  Because  of  the  quadratic 
dependency  of  local  ohmic  power  loss  on  current  density,  the  heat 
generation  in  different  cell  regions  will  grow  with  temperature.  As  a 
result  we  expect  a  higher  temperature  gradient  occurring  in  real  Li- 


ion  cells  compared  to  the  results  of  thermal  and  electric  cell-models 
that  disregard  the  influence  of  electric  and  thermal  interactions  in 
spatially  extended  systems. 

2.3.  Gradient  in  SoC 

In  the  following  it  is  demonstrated  that,  beside  a  current  density 
distribution,  a  SoC-difference  inside  a  LiFeP04-cell  is  expected  to 
occur  as  a  result  of  temperature  gradients  even  after  SoC-neutral 
current  load  profiles.  In  this  article  two  effects  will  be  discussed  as 
possible  reasons  for  SoC  inhomogeneities: 

1 )  The  varying  ratio  of  charge/discharge  impedance. 

2)  The  temperature  dependency  of  OCV. 

2.3.1.  SoC-differences  caused  by  varying  ratio  of  charge  I  discharge 
impedance 

As  we  can  see  in  Fig.  1,  the  pulse  impedance  varies  with  tem¬ 
perature  and  differs  for  charge  and  discharge  loads.  Fig.  2  shows 
that  additionally  the  ratio  of  discharge  to  charge  pulse  impedance 
varies  for  the  investigated  LiFeP04  cells  against  temperature.  In 
Fig.  2  the  ratio  of  the  10  s-discharge  pulse  impedance  to  the  1  Os- 
charge  pulse  impedance,  the  ZP(10s)-ratio,  is  plotted  against  SoC 
at  three  different  temperatures  (10,  23  and  40  °C). 

For  temperatures  between  23  °C  and  40  °C  the  ratio  of  discharge 
to  charge  pulse  impedance  decreases  with  increasing  temperature. 
In  contrast,  for  lower  temperatures  (23-1 0°C)  the  ratio  increases 
with  increasing  temperature  (except  for  25%  SoC).  These  qualita¬ 
tive  curves  were  found  for  the  ratio  of  the  30  s  pulse  impedance 
as  well.  An  explicit  electrochemical  reason  for  this  behavior  is  not 
yet  found  as  the  ratio  is  influenced  by  numerous  causes  and  as  the 
ratio’s  maxima  vs.  temperature  seems  to  be  the  result  of  at  least  two 
counteracting  thermodynamical  effects.  Flowever,  the  discussion 
of  the  experiments  will  show  a  good  agreement  of  the  measured 
ratios  for  the  observed  temperature  range  between  23  and  40  °C 
with  the  cells’  behavior  in  the  experimental  results. 

In  the  case  of  a  cycled  cell  with  temperature  gradient  inside, 
different  cell  regions  comprise  different  magnitudes  of  discharge 
to  charge  impedance  unbalances.  As  a  consequence,  current  den¬ 
sity  distributions  inside  the  jelly  roll  during  charge  and  discharge 
will  differ.  For  a  long-lasting  current  cycle  with  several  charge  and 
discharge  periods  these  unbalances  in  current  distribution  can  be 
considered  as  a  time-averaged  superposed  current  flow  from  one 
cell  region  to  another.  In  the  case  of  a  cell  operated  with  temper¬ 
ature  gradient  between  23  °C  and  40  °C  the  colder  cell  regions  of 
the  LiFeP04-cell  possess  larger  relative  discharge  impedance  and 
thus  will  be  less  discharged  during  discharge  mode  than  charged 
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Table  1 

Thermal  and  electrical  cell  specifications  of  the  investigated  Li-ion  cell. 
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Cathode  material 

LiFeP04 

Can  material 

Aluminum 

Anode  material 

Graphite 

Jelly  roll  length 

100  mm 

Nominal  voltage 

3.3  V 

Jelly  roll  density 

2680  kg  m-3 

Nominal  capacity 

4.4  Ah 

Jelly  roll  heat  capacity 

950JkgI<_1 

Geometrical  cell  type 

Cylindrical 

Jelly  roll  thermal  conductivity  (axial,  tangential) 

0.4  WmK-1 

Cell  housing  diameter 

32  mm 

Jelly  roll  thermal  conductivity  (radial) 

76WmI<-1 

Cell  housing  length 

113mm 

during  charge  mode.  Consequently,  even  during  SoC-neutral 
cycles,  the  local  SoC-values  inside  the  LiFeP04-cell  diverge.  This 
effect  (charge  drift  effect  1,  CDE  1)  is  limited  by  the  grow¬ 
ing  OCV  difference  caused  by  the  non-uniform  SoC  distribution 
in  the  cell.  As  a  result,  the  SoC  drift  should  converge  to 
a  limited  value  depending  on  the  charge  and  discharge  cur¬ 
rent  rates  of  the  load  profile.  At  other  temperature  ranges 
of  increasing  ZP-ratio  with  increasing  temperature  (below 
23  °C),  this  effect  results  in  a  SoC-difference  in  reversed  direc¬ 
tion. 

2.3.2.  SoC-differences  caused  by  the  temperature  dependency  of 
OCV 

The  second  effect  discussed  in  this  article  as  a  reason  for 
drifting  SoC  is  the  dependency  of  open  circuit  potential  on 
temperature,  which  can  be  seen  in  Fig.  3  for  the  investigated 
LiFeP04-cell.  Here,  voltage  profiles  were  monitored  for  cells  under 
open-circuit  conditions  during  variation  of  the  ambient  temper¬ 
ature  between  -20  and  40  °C.  These  measurements  were  done 
at  different  SoCs  reached  after  charge,  or  after  discharge  periods. 
Comparing  the  cells’  surface  temperature  profile  to  the  voltage 
profile,  the  linear  temperature  coefficient  of  OCV  (dD0/dT)  was 
determined.  For  50%  SoC,  the  positive  temperature  coefficient 
rises  up  to  0.1 2  mV  K-1.  Similar  values  for  LiFeP04 /graphite- 
cells  can  be  found  in  literature  [8].  The  shape  and  the  maxima 
of  the  dHo/dl  curve  against  SoC  correspond  with  the  entropy 
change  measurements  of  Viswanathan  et  al.  [9].  As  the  LiFeP04- 
electrode  seems  to  exhibit  only  small  changes  in  entropy  with 
varying  SoC,  the  graphite  electrode  can  be  identified  as  the  main 
influence  on  the  dD0/dT  characteristics.  According  to  Yazami 
and  Reynier  [10],  entropy  change-peaks  against  SoC  can  be 
explained  by  phase-transitions  during  lithium-ion  intercalation 
in  graphite  at  explicit  concentrations  of  lithium  atoms  in  Li*C6 
[10,11]. 
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Fig.  3.  Linear  temperature  coefficient  of  OCV  vs.  SoC  of  the  investigated  LiFeP04-cell. 


Transferred  to  the  given  problem  of  an  inhomogenously  tem¬ 
pered  cell,  warmer  cell  elements  reach  higher  OCV-values  than  the 
colder  elements  of  the  cell  assuming  the  same  SoC.  In  consequence 
of  this  potential  difference  between  the  cell  regions,  a  superposed 
transfer  current  occurs  discharging  the  warmer  cell  elements  with 
higher  OCV  and  charging  the  colder  elements.  The  current  flow 
remains  until  the  open  circuit  potential  difference  is  equalized.  The 
integrated  transfer  current  equates  the  fraction  of  SoC  drift  caused 
by  temperature  dependent  OCV.  Thus,  the  effect  of  growing  OCV 
against  temperature  (CDE  2)  also  leads  to  a  gradient  of  SoC  across 
the  elements  from  the  colder  and  higher  charged  elements  to  the 
warmer  and  lower  charged  elements. 

2.3.3.  Behavior  of  SoC-difference  during  relaxation  periods 

In  order  to  interpret  the  occurrence  and  decay  of  SoC-differences 
during  relaxation  periods  the  characteristics  of  the  open  circuit 
voltage  has  to  be  known.  In  the  first  case  we  assume  a  single-valued 
OCV-characteristic  against  SoC  (Fig.  4,  left).  This  means,  that  the 
OCV  of  a  cell  would  be  identical  after  charge  and  discharge  loads 
for  the  same  SoC  and  temperature.  Thus,  it  would  not  be  dependent 
on  the  load  history.  After  a  high  current  load  period  the  arisen  SoC 
drift  should  be  eliminated  directly  after  cycling  (CDE  1)  or  during 
temperature  equalization  (CDE  2),  respectively,  by  transfer  cur¬ 
rents  between  the  elements  due  to  potential  differences  (AOCV 
in  Fig.  4).  After  longer  rest  periods  the  SoC  difference  would  com¬ 
pletely  vanish.  This  is  schematically  depicted  in  Fig.  4  (left).  The 
SoCs  of  two  elements  in  a  cell— different  in  SoC  after  cycling  (A 
and  B)— assimilate  in  potential  and  thus  run  along  the  OCV-curve 
towards  each  other  to  the  same  threshold  (C).  No  difference  in  SoC 
would  persist.  This  assumption  is  a  good  approximation  for  elec¬ 
trode  chemistries  which  show  an  OCV-hysteresis  (see  below)  that 
is  small  compared  to  the  gradient  dD0/dSoC. 

Li-ion  cells  with  LiFeP04-cathodes  comprise  a  distinct  OCV  hys¬ 
teresis.  This  effect  is  characterized  in  several  publications  (e.g.  [12]) 
and  can  be  seen  as  the  coexistence  of  two  OCV-characteristics,  after 
charge  and  after  discharge  loads  in  consequence  of  varying  Li-ion 
concentrations  on  the  electrode  surfaces.  For  the  given  scenario  - 
relaxation  period  right  after  high  power  cycling  -  a  transfer  cur¬ 
rent  would  occur  similar  to  the  case  of  a  unique  OCV-characteristic 
degrading  the  SoC-difference.  If  a  higher  charge  throughput  ( ASoC 
>10%)  had  been  operated  right  before  the  end  of  cycling,  all  ele¬ 
ments  would  be  located  in  the  OCV-SoC  diagram  close  to  the  upper 
border  of  the  hysteresis  area  (A  and  B  in  Fig.  4  right). 

During  relaxation  for  the  two  observed  elements  in  Fig.  4  (right), 
the  element  with  lower  SoC  (A)  is  further  charged  by  the  transfer 
current  in  consequence  of  the  voltage  difference.  As  a  result  the 
element  runs  upwards  the  upper  OCV  curve  (C).  The  element  with 
higher  SoC  (B)  sees  a  change  in  the  current  direction  at  the  end 
of  cycling  and  gets  discharged  during  relaxation.  On  the  OCV-SoC 
diagram  the  point  of  element  B  moves  in  the  direction  of  the  lower 
OCV  curve,  the  discharge  OCV-curve  (C").  The  transfer  current  per¬ 
sists  until  all  elements  have  reached  the  same  OCV.  Considering 
the  non  unique  SoC-OCV  characteristic  a  SoC  difference  remains 
even  after  complete  compensation  of  the  element  temperatures 
(ASoCoo).  This  hysteresis  effect  influences  the  charge  balancing 
after  discharge  pulses  in  an  analog  way. 
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Fig.  4.  Behavior  of  the  local  OCV  of  a  high  charged  (A)  and  a  low  charged  (B)  jelly  roll  element  during  relaxation  in  a  schematic  OCV-SoC-diagram,  left:  considering  a  unique 
OCV-SoC  characteristic,  resulting  in  a  single  SoC-equalization  (C);  right:  considering  a  OCV  hysteresis  in  dependency  on  current  direction.  The  OCV  of  B  runs  during  relaxation 
to  point  C”.  The  original  SoC-difference  ( ASoCo)  is  reduced  to  a  constant  SoC-difference  ( ASoCoo).  For  further  explanation  see  Section  2.3.3. 


3.  Simulation  model 

It  is  well-known,  that  the  interdependencies  between  thermal 
and  electrochemical  behavior  of  a  battery-cell  make  a  permanent 
co-simulation  of  an  electric  and  a  thermal  model  mandatory  to 
predict  the  cell’s  behavior  precisely.  This  was  considered  by  Gu 
and  Wang  [13],  who  coupled  the  volume-averaged  thermal  energy 
equation  with  a  first-principle  electrochemical  model.  This  model 
uses  the  heat  generation  and  the  temperature-dependent  phys- 
iochemical  properties  of  transport,  kinetics  and  mass-transfer  as 
interfaces  between  thermal  and  electric  model.  Later  Srinivasan 
and  Wang  [14]  extended  a  similar  model  with  a  thermal  2D-model 
including  the  local  heat  generation  method. 

The  following  thermal-electric  simulation  model  of  the  cylin¬ 
drical  LiFeP04-cell  was  used  to  reproduce  the  experimental 
part  and  on  the  other  hand  in  the  later  discussion  to  quantify 
thermal-electric  inhomogeneities  inside  a  single  LiFeP04-cell  dur¬ 
ing  operation.  Therefore,  local  thermal  and  electrical  distributions 
inside  the  cell  had  to  be  calculable.  Thus,  a  spatially  resolved  three 
dimensional  thermal  model  in  combination  with  an  impedance 
based  spatially  resolved  three  dimensional  electric  model  of  the 
cell’s  jelly  roll  is  used  for  simulation.  The  time-discrete  model  was 
implemented  using  the  simulation  software  Matlab  Simulink  [15]. 
Detailed  information  on  the  two  submodels  is  given  in  the  follow¬ 
ing,  and  was  also  reported  previously  in  [16]. 


tion,  4  in  tangential  direction).  Thus,  the  discretization  divides  the 
jelly  roll  in  140  discret  elements  of  the  same  volume.  The  temper¬ 
ature  of  each  volume  element  (VE)  is  calculated  by  thermal  energy 
conservation  considering  reversible  and  irreversible  heat  sources 
of  each  VE  and  interacting  heat  fluxes  between  the  adjoining  VEs. 
The  explicit  numeric  solution  of  the  differential  equation  for  one 
VE  in  order  to  calculate  the  temperatures  for  the  next  time  step  out 
of  the  given  temperature  distribution  is  given  in  Eq.  (3). 

TAX  +  At)  .  Tilt)  +  ^  ^  +  Oirr.E  +  <w)  At 

(3) 

Subscript  E  (for  element)  indicates  a  VE  specific  value  (tempera¬ 
ture  or  heat  source).  Subscript  n  indicates  the  number  of  adjoining 
VEs,  Ti  is  the  temperature  of  the  ith  adjoining  VE  and  ftthiEl-  is 
the  thermal  resistance  from  VE  E  to  the  ith  adjoining  VE.  Rt h  is 
a  geometry  dependent  value  influenced  by  thermal  conductivity. 
The  thermal  conductivity  and  the  specific  heat  capacity  of  the 
homogenous  jelly  roll  are  assumed  to  be  constant  in  the  investi¬ 
gated  temperature  range  and  independent  of  SoC. 

Irreversible  and  reversible  heat  generation  for  each  control 
volume  is  mathematically  implemented  in  the  cell  simulations 
according  to  Wu  et  al.  [18]: 


3.1.  Thermal  model 


The  correlation  of  heat  flux  and  temperature  development 
inside  the  jelly  roll  depends  on  the  phenomenon  of  heat  conduc¬ 
tivity  and  is  described  by  Fourier’s  law  of  heat  conduction  [17]. 
Considering  the  cell’s  geometry,  it  is  applied  for  cylindrical  coordi¬ 
nates: 


Av2r+^  =  |  with  v2t 
pc  pc  St 


1  8_  (  ST\  \  S^T 
r  Sr  l  r  Sr  J  +  r 2  Sep 2  +  Sz2 

(2) 


Abbreviations  and  parameter  definitions  of  all  equations  are 
given  in  Table  Al. 

The  jelly  roll,  consisting  of  a  spirally  wound  electrode  stack,  is 
macroscopically  modeled  as  a  homogenous  hollow  cylinder  with 
anisotropic  heat  conduction  properties  -  exhibiting  a  low  radial  and 
a  high  axial  and  tangential  thermal  conductivity.  The  differential 
equation  of  heat  conduction  is  discretizised  with  the  finite-volume- 
method  (FVM)  that  is  applied  for  three  dimensions  in  cylindrical 
coordinates  (7  volume  elements  in  radial  direction,  5  in  axial  direc¬ 


n 

Qirr  =  ^.E  X  lf  E  (4) 

1  =  1 

Qrev  =  TE^x/E  (5) 

Rt  indicates  the  ith  ohmic  resistance  of  a  user-defined  equivalent 
circuit  model  (ECM)  with  n  ohmic  resistances  inside  the  element  E. 
/j)E  indicates  the  current  flow  through  Ri  E.  U0,E  is  the  open  circuit 
voltage  of  element  E  and  an  approximation  for  the  thermodynamic 
equilibrium  voltage  that  can  be  seen  as  the  threshold  of  the  OCV 
for  infinitely  long  resting  periods.  /E  is  the  current  flow  through  the 
serial  elements  of  the  used  local  ECM  in  element  E.  The  heat  gener¬ 
ation  effect  by  heat  of  mixing  and  parasitic  chemical  reactions,  as 
discussed  in  [1 8,1 9],  do  play  a  minor  role  in  overall  heat  generation 
during  operation  and  are  neglected. 

The  cell  peripherals  -  anode  and  cathode  terminals  and  housing 
-  were  modeled  as  a  network  of  thermal  capacities  and  resistances. 
The  thermal  resistances  of  free  convection  on  the  cell  housings’ 
surface  are  considered  by  a  global  heat  transfer  coefficient  a  on  cells 
surface.  The  values  of  a  were  determined  in  previous  experiments. 
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ume  element  with  its  local  OCV  (Uo.e)  and  the  ohmic  resistances  of  the  anode  and 
cathode  current  collector  to  the  adjoining  volume  elements  (in  axial  and  tangential 
direction). 

3.2.  Electrical  model 

For  the  reproduction  of  the  electric  cell  behavior  an  equiva¬ 
lent  circuit  model  (ECM)  is  employed.  Buller  [20]  gives  a  review 
on  differently  complex  ECMs  used  in  Li-ion  cell  simulations.  The 
model  chosen  for  the  described  investigation  consists  of  a  serial 
resistance  in  combination  with  two  RC-circuits  in  series.  All  ele¬ 
ments  are  parameterized  as  a  function  of  SoC,  temperature  and 
current  direction.  The  values  were  assumed  to  be  independent 
from  current  in  the  investigated  current  range.  As  depicted  in  Fig.  5 
one  RC-circuit  describes  charge  transfer  resistance  (RCt)  and  dou¬ 
ble  layer  capacity  (Cdl),  the  other  describes  the  impedance  of  the 
solid  electrolyte  interface  (RSei  and  Cs ei)  according  to  [21].  The 
serial  resistance  (Rs)  is  mainly  influenced  by  the  limited  conduc¬ 
tance  of  contacts,  active  materials  and  electrolyte.  Cell  inductivity 
and  diffusion  overpotentials  caused  by  concentration  gradients  - 
often  modeled  by  a  Warburg  impedance  -  were  neglected.  The 
relatively  simple  ECM  structure  is  a  compromise  between  accu¬ 
racy  and  computing  time,  as  the  spatialization  of  the  electric  cell 
model  requires  multiple  solving  of  the  ECM  per  time  step  -  once 
for  each  VE.  A  method  for  implementing  various  components  of 
ECMs  in  Matlab  Simulink  for  simulation  is  also  given  by  Buller 
[20]. 

As  shown  in  Fig.  5  the  described  ECM-model  of  the  complete 
cell  was  spatialized  to  an  electric  Finite-Network-Method  (FNM- 
)  model,  in  which  each  volume  element  of  the  jelly  roll  (index  E 
in  Fig.  5)  exhibits  its  own  local  anode  potential,  cathode  potential, 
local  impedance,  local  OCV,  local  heat  generation,  and  local  SoC. 
For  a  distinct  thermal-electric  co-simulation  the  FNM-model  was 
discretized  according  to  the  thermal  model.  Thus,  the  nodes  on  the 
cathode  current  collector  of  each  VE  are  connected  to  the  nodes 
of  the  adjoining  VEs  by  pure  ohmic  resistances  in  tangential  (Rta,c) 
and  axial  direction  (Rax,c)  calculated  by  material  and  geometrical 
properties.  In  analogy,  this  was  implemented  for  the  anode  nodes 
(Rta  A  and  Rax>A).  The  axially  outer  elements  of  each  radial  layer  of 
VE  are  connected  to  the  accordant  terminal  by  ohmic  resistances  of 
the  terminal  current  collectors.  These  resistances  were  estimated 
by  geometry  and  material  properties.  The  electric  resistances  of 
anode  and  cathode  current  collector,  as  well  as  the  terminal  cur¬ 
rent  collectors,  are  at  least  one  order  of  magnitude  smaller  than 
the  serial  resistance  of  the  local  ECMs.  Thus  the  estimated  ohmic 
resistances  have  little  impact  on  the  simulation  results. 

Additionally  for  every  element  a  linear  hysteresis  model  for 
determining  OCV  as  a  function  of  SoC,  temperature,  and  load 
history  is  integrated  -  a  simplification  of  the  OCV-model  of 


Thele  et  al.  [22]. 

Uo,E  =  Y  ^/0,ch(SoCE)  +  ±^(rE  -  25  °C) 

+(1  -  Y)  ^/0,dch(SoCE)+  -25°C)j  (6) 

The  OCV-curves  of  U0  ch  and  Uo  dch  are  implemented  in  3D  look 
up  tables  in  dependence  of  SoC  and  temperature.  The  coefficient 
7  describes  the  linear  hysteresis  model.  It  reaches  its  maximum 
value  of  1  after  long  charging  periods  (>5%  of  SoC)  and  minimum  0 
after  long  discharging  periods  (>5%  of  SoC).  Otherwise  it  decreases 
linearly  with  discharge  and  increases  linearly  with  charge  current 
throughput  between  0  and  1  with  a  gradient  of  l/(0.05  x  Cn0m). 


4.  Experimental 


Measuring  temperature  distributions,  local  current  distribu¬ 
tions  or  SoC-gradients  inside  a  jelly  roll  is  not  possible  in  an 
experimental  set-up  without  disturbing  cell  configuration  and  per¬ 
formance.  As  an  indirect  evidence  of  local  electric  differences, 
thermal-electric  interactions  have  been  validated  by  observing 
three  cells  on  different  temperature  levels  electrically  connected 
in  parallel.  This  test  was  also  modeled  in  thermal-electric  cell 
simulations  for  validation  of  the  applied  model.  Object  of  the  inves¬ 
tigations  were  cylindrical  LiFeP04  high  power  cells  (specification  in 
Table  1). 

Firstly,  a  cylindrical  Li-ion  cell  with  a  continuous  temperature 
and  current  density  distribution  can  simply  be  modeled  as  three 
radially  arranged  volume  elements  electrically  connected  in  paral¬ 
lel.  The  three  volume  elements  exhibiting  different  temperatures 
represent  the  inner,  middle  and  outer  region  of  the  cell.  Further¬ 
more,  each  of  the  three  elements  is  considered  in  the  experiments 
as  single  cell  (Fig.  6,  middle).  In  order  to  rebuild  temperature  gra¬ 
dients  between  the  elements,  the  three  cells  connected  in  parallel 
were  operated  at  different  temperatures.  The  thermal  differences 
of  the  three  cells  during  operation  were  induced  by  different  heat 
dissipation  paths  on  the  cells’  surfaces.  Cell  1  was  radially  thermally 
insulated  by  a  20  mm  thick  plastic  foam,  cell  2  was  encased  by  an 
insulation  of  plastic  foam  as  well,  but  only  for  half  of  its  length, 
and  cell  3  was  exposed  to  free  convection  in  contact  with  5°C 
cooled  ambient  air.  The  current  flow  through  every  cell  was  mea¬ 
sured  independently  by  shunts  connected  in  series  to  the  cells.  The 
three  electric  paths  were  routed  through  a  contactor  to  electrically 
separate  the  cells  after  operating.  Different  sensors  monitored  volt¬ 
age  (absolute  accuracy:  ±5  mV),  current  (accuracy:  ±0.05  A)  and 
housing  temperature  (accuracy:  ±1  K)  of  every  cell. 

Fig.  6  (right)  shows  the  schematic  diagram  of  the  electric  circuit 
that  was  connected  to  a  battery  test  bench.  All  experiments  were 
executed  in  a  climatic  chamber. 

Previous  to  the  main  experiments  (ME  1 ),  the  three  most  similar 
cells  were  chosen  from  a  set  of  six  cells  by  comparing  their  actual 
capacity  and  their  internal  impedance.  The  relative  differences  of 
measured  capacity  between  the  three  chosen  cells  were  below  0.5%, 
the  differences  of  the  1  OC-pulse  impedance  below  1  %.  All  measured 
SoCs  in  this  article  are  the  ratio  of  the  actual  usable  charge  in  cell 
for  lC-discharge  current  (/n0m)  at  25  °C  and  a  cutoff  voltage  of  2  V 
to  the  nominal  capacity. 


SoC  = 


f 

Jo 


t(U=2V) 


fnom  dt 


T=25°C 


(7) 


In  order  to  start  the  experiment  with  all  three  cells  at  equal 
states,  the  cells  in  the  setup  were  first  cooled  down  by  5  °C  ambient 
temperature  and  then  completely  discharged  in  parallel  (cutoff- 
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(represents  inner  cell  region) 
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(represents  middle  cell  region) 
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Single  cell  Model  of  three  Experimental 
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•  Thermocouple  type  K 
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Fig.  6.  Experimental  approach  of  three  cells  at  different  temperatures  electrically  connected  in  parallel;  left:  transformation  of  a  single  cell  with  continuous  temperature 
distribution  to  a  model  with  three  radial  arranged  volume  elements;  middle:  representation  of  the  three  VEs  by  the  three  cells;  right:  schematic  experimental  setup  with 
contactor  and  shunts  connected  to  the  battery  test  bench. 


voltage  2  V).  Hereupon  they  were  charged  with  low  current  (2/3C) 
to  62%  SoC  and  then  rested  for  30  m  at  5°C.  For  an  initial  warm¬ 
up  the  cells  were  first  heated  by  5C  pulse  cycles  (5C  for  each  cell: 
3  x  5C  =  66 A)  for  15m  (90s  discharge-90 s  charge).  Then  the  cur¬ 
rent  amplitude  was  increased  to  3  x  8C  (105.6  A)  and  the  cells  were 
operated  for  another  2  h.  The  cells  reached  a  thermal  steady  state 
during  this  period.  At  the  end  of  the  SoC-neutral  cycle,  the  cells 
were  immediately  separated  by  opening  the  contactors.  After  tem¬ 
pering  the  cells  to  room  temperature,  their  SoCs  were  separately 
determined  by  lC-discharging  to  a  cutoff- voltage  of  2  V. 

In  a  further  experiment  (ME  2)  ME  1  was  exactly  repeated 
without  separating  the  cells  after  cycling.  Instead  transfer  cur¬ 
rents  during  relaxation  were  enabled  and  measured  for  10  h  at 
23  °C  ambient  temperature.  This  time  period  was  significantly 
longer  than  the  cooling-down-period  of  all  three  cells.  Again,  SoC- 
determinations  of  all  cells  were  performed  after  separating  the  cells 
analog  to  the  previous  tests. 

For  the  simulative  reproduction  of  the  described  experiments 
the  heat  transmission  coefficients  of  the  cell  surfaces  to  environ¬ 


ment  were  measured  by  previous  experiments  as  the  undefined  air 
flowing  conditions  in  the  climate  chamber  prevent  an  exact  calcu¬ 
lation  of  these  values.  Therefore  the  three  cells  were  operated  in  the 
experimental  setup  of  ME  1  with  SoC-neutral  current  pulse  cycles 
(90  s  discharge-90  s  charge)  and  changing  current  amplitude,  each 
time  until  a  thermal  steady  state  was  reached.  The  cells’  power  loss 
was  estimated  according  to  Onda  et  al.  [23].  Considering  the  dif¬ 
ference  of  cell  surface  to  ambient  temperature  the  heat  transfer 
coefficients  were  calculated  to  ot\  =15.5,  a2  =  23.9,  a3  =40.1  using 
Eq.(8). 


(7s,avg  -  Tamb)As 

5.  Results 

The  temperature  gradient  between  the  three  cells  of  maximal 
20 1<  (fully-isolated  to  non-isolated)  can  be  seen  in  Fig.  7.  Herein 
the  measured  and  predicted  temperature  profiles  of  the  three  cells 


Fig.  7.  Temperature  development  during  the  pulse  cycle  of  ME1  for  the  three  cells  comparing  simulated  (dashed  lines)  and  measured  values  (solid  lines);  blue:  non  isolated 
cell,  green:  semi  isolated  cell,  red:  isolated  cell.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


Table  2 

Comparison  of  measured  and  simulated  average  current  values  for  the  three  cells  during  the  last  charge  and  discharge  current  pulse  at  the  end  of  ME  1  (isolated  cell:  1, 
semi-isolated  cell:  2,  non-isolated  cell:  3). 


Average  current  of  last  charge  pulse  /Ch,avg  (A) 

Average  current  of  last  discharge  pulse  /dchjavg  (A) 

Measured 

Simulated 

Measured 

Simulated 

Cell  1 

40.4 

40.1 

-40.5 

-40.1 

Cell  2 

34.6 

34.6 

-34.6 

-34.5 

Cell  3 

30.1 

30.1 

-30.3 

-31.0 
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Fig.  8.  Current  development  during  the  thermal  steady  of  the  last  cycle  pulses  of  ME  1  for  the  three  cells  simulated  (dashed  lines)  and  measured  values  (solid  lines);  blue: 
non  isolated  cell,  green:  semi  isolated  cell,  red:  isolated  cell.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the 
article.) 


on  their  surfaces  are  plotted  against  time.  A  good  agreement  of 
simulation  and  experimental  curves  can  be  seen.  The  difference  of 
maximal  2 1<  (occurring  for  the  fully  isolated  cell)  can  be  explained 
by  the  variance  of  the  isolation  efficiency  between  previous  and 
main-experiment. 

Caused  by  growing  temperature  difference  the  fraction  of  cur¬ 
rent  flow  through  the  isolated  cell  in  relation  to  the  entire  current  is 
expected  to  increase,  until  thermal  steady  state  is  reached.  The  good 
agreement  of  simulation  and  experiment  can  be  seen  on  the  cur¬ 
rent  diagram  of  the  three  cells  for  the  last  pulses  of  cycling,  where 
thermal  steady  state  is  already  reached  (Fig.  8).  The  average  cur¬ 
rents  for  the  three  cells  during  the  last  charge  (/Ch,avg) and  discharge 
pulse  (/dch.avg) are  compared  in  Table  2.  The  parasitic  resistances  of 
the  connecting  elements  and  shunts  influencing  the  current  dis¬ 
tribution  were  measured  to  2.15  ±0.04  m£2  and  considered  in  the 
simulation. 

For  the  analyzed  test  the  fully  isolated  cell  had  -  as  expected  - 
the  lowest  SoC  right  after  disconnecting  (SoCend0)  of  the  three  cells 
while  the  non-isolated  cell  was  the  highest  charged.  The  values  can 
be  seen  in  Table  3.  The  impact  of  OCV  hysteresis  of  the  LiFeP04- 
cell  on  SoC-behavior  during  relaxation  periods  can  be  illuminated 
by  comparing  the  results  of  ME1  and  ME2.  For  the  experiment 
without  immediate  electric  separation  (ME  2)  we  can  assume,  that 
the  cells’  SoCs  directly  after  operation  were  similar  to  the  values 
measured  for  ME1  (SoCend  0).  The  SoC  of  the  cells  after  the  resting 
period  of  lOh  (SoCendoo)  were  determined  by  1C  discharge  analog 
to  ME1.  The  difference  between  SoCend  0  and  SoCend  oo  for  each  cell 
represents  its  SoC-compensation  during  relaxation  (ASoCtransfer). 
To  check  the  values  of  SoC-compensation  the  cell  currents  were 
recorded  and  integrated  during  relaxation  period.  Both  methods  of 
determining  ASoCtransfer  led  to  similar  results. 

It  can  be  seen,  that  during  longer  rest-periods  after  high  current 
loads  the  SoCs  of  Li-ion  cells  with  hysteresis  behavior  connected  in 
parallel  does  not  equalize  completely.  For  the  investigated  oper¬ 
ation  only  about  one  third  of  the  SoC-difference  was  equalized 
during  operation  both,  in  the  experiment  and  in  simulation.  The 
cells  persist  stably  at  differing  SoC-values. 

6.  Discussion 

The  experimental  results  in  accordance  with  the  simula¬ 
tions  detect  thermal-electric  inhomogeneities  between  a  set 
of  parallel  connected  LiFeP04-cells  of  differing  temperature. 
These  inhomogeneities  have  to  be  expected  in  non-uniformly 
cooled  battery-packs  containing  cells  connected  in  parallel.  The 
differences  in  temperature,  current  density  and  SoC  lead  to  inho¬ 
mogeneous  cell  aging  inside  the  battery-pack  influencing  the 


overall  battery  performance  in  the  course  of  time.  Although 
the  results  may  easily  be  transferred  to  multi-cell  setups,  the 
focus  of  the  following  discussion  is  the  quantification  of  the 
thermal-electric  interactions  that  are  expected  inside  a  single  cell. 

A  LiFeP04-cell  -  exposed  to  forced  convection  -  was  simu¬ 
lated  using  the  spatialized  thermal/electric  cell-model.  Therefore 
the  ambient  temperature  of  5°C  and  a  SoC-neutral  pulse-cycle 
with  a  magnitude  of  70  A  were  chosen  (20  s  discharge  pulses-20  s 
charge  pulses).  The  starting  SoC  was  set  to  50%.  The  heat  transfer 
coefficient  for  an  incident  flow  (lateral  to  the  middle  axis)  with  a 
velocity  of  15 ms-1  and  a  cylindrical  geometry  was  calculated  to 
1 1 7.9  W m2  K-1  on  cell  surface  according  to  [24]. 

Fig.  9  show  snapshots  of  the  respective  cell  state  distributions 
at  the  end  of  the  last  charge  pulse  (after  6530  s)  in  a  thermal 
steady  state.  All  three  depictions  contain  two  sectional  views,  on 
the  left  side  in  the  radial-axial  plane  and  on  the  right  side  in 
radial-tangential  plane  (in  the  axial  middle  of  the  cell). 

The  simulated  temperature  distribution  in  the  jelly  roll  can  be 
seen  in  Fig.  9a.  The  maximum  temperature  of  jelly  roll  is  40.9  °C 
at  the  inner  part  of  the  jelly  roll.  The  negative  temperature  gradi¬ 
ent  increases  to  the  outer  regions  of  the  jelly  roll.  The  minimum 
temperature  in  the  jelly  roll  for  the  given  scenario  is  19.5  °C.  Fig.  9b 
shows  the  distribution  of  current  flow  (between  anode  and  cath¬ 
ode)  per  volume  (current  density)  with  a  similar  shape  as  the 
temperature  distribution.  It  can  be  seen  that  the  current  density 
in  the  inner  and  outer  parts  of  the  jelly  roll  differ  by  a  factor  of 
approx.  1.3,  where  higher  currents  are  recorded  for  the  inner  cell 
region.  The  variation  of  local  SoC  is  monitored  in  Fig.  9c.  For  this 
point  in  time,  where  the  average  cell  SoC  is  equal  to  the  starting 
SoC,  the  maximal  SoC-difference  reaches  9%  ( ASoCend,o). 

Fig.  10  shows  the  SoC  profile  of  the  hottest  and  coldest  region 
inside  the  jelly  roll  vs.  time.  The  two  exemplary  local  SoCs  dis¬ 
perse  induced  by  rising  temperature  difference  and  do  not  equalize 
completely  after  end  of  cycling  according  to  the  influence  of  OCV- 
hysteresis  (ASoCendloo  =  ca-  5%).  According  to  the  simulation  a  cell 
that  is  stored  after  high  current  loads  would  possess  a  SoC-gradient 
persisting  during  rest  periods. 

In  summary  the  three  following  spatial  thermal-electric  effects 
are  expected  to  influence  the  aging  mechanism  and  thus  the  life¬ 
time  of  LiFeP04-cells  as  the  aging  at  different  jelly  roll  locations 
proceeds  under  different  aging  condition:  generally  the  impact  of 
temperature  on  calendar  life  predictions  of  Li-ion  cells  is  mostly 
considered  as  an  exponential  dependence  (Arrhenius  equation) 
[5,25,26].  Also,  an  exponential  aging  behavior  vs.  temperature  is 
reported  in  several  capacity  fading  investigations  at  cycle  life  con¬ 
ditions  [1,5].  For  LiFeP04/C-cells  the  effect  of  accelerated  cycle  life 
aging  on  elevated  temperatures  can  also  be  seen  in  [6].  Assuming 
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Fig.  9.  Simulated  temperature  distribution  (a),  current  density  distribution  (b)  and  SoC-distribution  (c)  inside  the  LiFeP04-cell  directly  at  the  end  of  the  single  cell  cycle 
(t=  6530s). 


an  Arrhenius  relation,  a  temperature  gradient  inside  a  cell  induces 
an  exponentially  rising  aging  behavior  of  the  inner  regions  of  a 
cylindrical  Li-ion  cell. 

The  current  distribution  caused  by  the  temperature  gradient  is 
expected  to  intensify  this  unequal  aging  behavior.  E.g.  Ning  et  al. 
assumes  that  for  LiCo02 -cells  capacity  fade  approx,  grows  with 
the  square  root  of  the  current  magnitude  [27].  Additionally  the 
increased  charge  throughput  in  warmer  cell  regions  is  expected 
to  induce  locally  accelerated  capacity  fading.  According  to  charge 
throughput  models,  which  are  state  of  the  art  for  empirical  cycle  life 


aging  models  [28],  aging  progress  is  in  direct  correlation  to  charge 
throughput.  Thus,  a  further  increased  gradient  of  aging  behavior  in 
the  interior  regions  of  the  cell  is  expected. 

Also  thermally  induced  SoC-distributions  can  influence  the 
gradient  of  aging  progresses.  From  other  cell  chemistries,  the 
SoC-dependence  of  cell  aging  is  well  known.  For  example  Amine 
et  al.  [29]  showed  a  significant  difference  in  power  fading  for  a 
LiNiCo02/C-cell  stored  at  40%  and  60%  SoC  at  several  tempera¬ 
tures.  Additionally,  faster  capacity  degradation  is  expected  to  occur. 
However,  for  LiFeP04/C-systems,  which  tend  to  significant  SoC- 


Table  3 

Comparison  of  cells’  SoCs  at  the  end  of  cycling  (SoCend,o)  and  after  10 h  resting  period  (SoCend,oo)  in  simulation  and  in  experiment.  SoC-compensation  during  relaxation 
( ASoQransfer)  is  calculated  firstly  by  the  difference  of  SoCend,oo  and  SoCend.o  and  secondly  by  current  integration  during  relaxation. 


SoCend,0 

SoCend>00 

ASoQransfer 

Measured 
(ME  1) 

Simulated 

Measured 
(ME  2) 

Simulated 

Calculated 

Current 

integrated 

Cell  1 

59.1% 

58.9% 

60.1% 

60.0% 

+1.0% 

+0.9% 

Cell  2 

63.1% 

63.3% 

62.6% 

62.7% 

-0.6% 

-0.4% 

Cell  3 

64.4% 

63.9% 

63.7% 

63.3% 

-0.7% 

-0.6% 

ASoC 

-5.3% 

-5.0% 

-3.6% 

-3.3% 

Fig.  10.  SoC  progress  of  the  hottest  (blue  line)  and  coldest  cell  region  (green  line)  inside  the  LiFeP04-cell  during  the  single  cell  cycle  with  a  persisting  SoC  difference 
(ASoCend.oo)-  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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Table  Al 

Abbreviation  and  parameter  definition  used  in  formulas,  text  and  figures. 


4777 


Zp(t) 

Pulse  impedance  [£2] 

Rs 

Serial  resistance  of  the  ECM  [£2] 

A  U 

Internal  voltage  drop  of  cell  [V] 

Rc  T 

Charge  transfer  resistance  of  the  ECM  [Q] 

Al 

Current  step  magnitude  [A] 

Cdl 

Double  layer  capacity  of  the  ECM  [F] 

Ho 

Open  circuit  voltage  [V] 

Rsei 

Resistance  of  the  solid  electrolyte  interface  in  the  ECM 

[ni 

Capacity  of  solid  electrolyte  interface  in  the  ECM  [F] 

T 

Temperature  [K] 

CsEI 

X 

Thermal  conductivity  [W  ml<-1  ] 

Rax,C 

Ohmic  resistance  of  cathode  current  collector  in  axial 
direction  between  two  CV  [Q,] 

P 

Density  [kg/m3] 

Rta,C 

Partial  resistance  of  cathode  c.c.  in  tan.  direction  [Q,] 

c 

Specific  heat  capacity  [J  kg  K-1  ] 

RaxA 

Partial  resistance  of  anode  c.c.  in  axial  direction  [£2] 

<7q 

Density  of  heat  source  [  W  m-3  ] 

R-ta.A 

Partial  resistance  of  anode  c.c.  in  tan.  direction  [£2] 

t 

Time  [s] 

Ho.ch 

OCV  after  charge  load  history  [V] 

r 

Radial  coordinate 

Ho.dch 

OCV  after  discharge  load  history  [V] 

z 

Axial  coordinate 

7 

Linear  OCV-hysteresis  variable  [-] 

(P 

Tangential  coordinate 

Cnom 

Ij iom 

Nominal  capacity  [Ah] 

Constant  discharge  current  of  4.4A  (1C)  [A] 

At 

Time  step  size  of  simulation  [s] 

a 

Heat  transfer  coefficient  [W(m2I<)_1  ] 

VE 

Volume  of  control  element  [m3] 

Pv.avg 

Average  cell  power  loss  [W] 

Rth 

Thermal  resistance  [K  W_1  ] 

Lamb 

Ambient  temperature  [K] 

Qirr 

Irreversible  heat  generation  [W] 

Ls.avg 

Average  cell  surface  temperature  [K] 

Qrev 

Reversible  heat  generation  [W] 

As 

Cell  surface  area  [m2] 

Ri 

ith  ohmic  resistance  of  an  user-defined  equivalent 
circuit  model  [£2] 

Lh.avg 

Average  cell  current  during  the  last  charge  pulse  (ME 

1)  [A] 

h,E 

Current  flow  through  the  ith  ohmic  resistance  in  the 

ECM  of  element  E  [A] 

fdch.avg 

Average  cell  current  during  the  last  discharge  pulse 
(ME  1)  [A] 

h 

Current  flow  through  the  serial  elements  in  the  ECM  of 
Element  E  [A] 

SoCend.O 

Cell’s  SoC  at  end  of  cycling 

Lrell 

Temperature  of  cell  [°C] 

SoCend,oo 

Cell’s  stable  SoC  after  resting  period 

Lell 

Cell  current  [A] 

ASoCtransfer 

Eliminated  SoC  difference  during  relax  period. 

gradients,  the  magnitude  of  SoC-influence  on  calendar  and  cycle  life 
aging  -  and  thus  on  inhomogeneous  aging  inside  a  single  cell  -  has 
to  be  investigated  in  future.  Furthermore,  different  electrochemical 
aging  effects  that  are  influenced  by  high  temperature,  high  cycling 
rates  and  high  and  low  SoCs  in  Li-ion  cells  are  discussed  by  Vetter 
et  al.  [3]. 

In  summary,  deduced  from  literature  the  three  observed 
inhomogeneities  (temperature,  current  density  and  SoC)  can,  in 
common,  lead  to  a  more  than  linearly  accelerated  aging  progress  for 
the  warmer  inner  regions  inside  a  Li-ion  cell  compared  to  the  colder 
outer  regions.  In  consequence,  a  different  and  faster  aging  behavior 
of  the  entire  cell  is  expected  than  for  a  cell  with  -  hypotheti¬ 
cally  -  absolutely  uniform  temperature  distribution  with  the  same 
volume  averaged  temperature.  Thus,  the  consideration  of  local 
thermal-electric  interactions  and  the  usage  of  spatially  resolved 
thermal-electric  cell  models  are  hence  important  for  aging  analy¬ 
sis  and  therefore  for  thermal  management  aspects  of  Li-ion  cells  in 
order  to  optimize  active  cooling  concepts  and  strategies. 

7.  Conclusion 

In  this  article  the  thermal-electric  interactions  caused  by  tem¬ 
perature  gradients  inside  a  LiFeP04  cell  are  treated.  In  addition  to 
temperature  and  current  gradients  a  divergence  of  the  local  state  of 
charge,  the  SoC  drift,  could  be  explained  and  verified  experimen¬ 
tally  even  for  charge  neutral  current  profiles.  The  following  two 
effects  could  be  identified  as  the  reasons  of  the  SoC-drift: 

1.  The  temperature  dependent  ratio  of  charge  to  discharge  pulse 

impedance  (CDE  1). 

2.  The  temperature  dependence  of  open  circuit  potential  (CDE  2). 

Electric  inhomogeneities  between  three  differently  cooled  and 
electrically  in  parallel  connected  LiFeP04-cells  were  detected  both, 
in  experiment  and  simulation.  In  the  experiments  the  maximum 
difference  in  SoC  of  5.3%  after  a  SoC-neutral  pulse  cycle  in  a  ther¬ 
mal  steady  state  occurred  between  two  cells  with  a  difference 
in  temperature  of  20  K.  The  results  of  simulation  agree  well  with 


temperature,  current  and  SoC  measurements.  Both  examinations 
show  that  for  the  given  experiment  even  after  long  resting  periods 
the  SoC-difference  between  the  cells  is  reduced  only  to  2/3  of  the 
original  value  (5.3  -*  3.7%).  An  SoC  difference  persists  due  to  the 
OCV-hysteresis,  which  is  significant  in  LiFeP04/C-cells. 

The  conclusions  for  cells  connected  in  parallel  were  additionally 
transferred  to  the  single  cell  behavior.  Spatially  resolved  ther¬ 
mal  and  electric  cell  co-simulations  of  a  LiFeP04  cell  loaded  with 
a  70  A  pulse  cycle  showed  a  current  density  gradient  across  the 
jelly  roll  (iinner/Wer:  ca.  1.3)  as  well  as  a  SoC-incline  of  9.5% 
from  inner  to  outer  jelly  roll  elements  caused  by  a  temperature 
spread  from  19.5  °C  (outer  elements)  to  40.9  °C  (inner  elements). 
Thus  temperature,  current  and  SoC  distribution  indicates  more 
harsh  aging  conditions  for  the  warmer  parts  of  a  Li-ion  cell  as 
well  as  an  accelerated  aging  of  the  complete  cell  in  compar¬ 
ison  to  cells  with  homogenous  temperature  distributions.  The 
importance  to  study  the  thermal-electrical  inhomogeneities  for 
accurate  aging  examinations  can  be  deduced  from  the  literature. 
Investigations  and  experiments  considering  aging  inhomogeneities 
have  to  be  done  both  on  single  cell  and  on  battery  pack 
level. 

Appendix  A. 

See  Table  Al. 
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